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of the electroweak gauge bosons self-energies for arbitrary momentum transfer 
and for different internal quark masses. Compact formulae for both the real 
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is discussed. 



*NSERC Fellow. 



1 



1. Introduction 



The Standard Model of the electroweak and strong interactions has achieved a tremendous 
success in describing all experimental data within the range of energies available today. In 
particular, it has been tested to the level of its quantum corrections in the high-precision LEP 
and SLC experiments 0. To allow for such precision tests, a considerable amount of theoretical 
work has been done in the last few years to calculate the relevant radiative corrections [[J. 
In fact, in view of the remarkably high accuracy of the experimental data, some radiative 
corrections were required beyond the one-loop approximation. 

Except for the QCD corrections to light quark pair production in e + e~ annihilation, which are 
known up to 0(a 2 s ) |J since more than a decade and which have been recently calculated to 
0{pP s ) PJ, these two-loop corrections were mainly concerned with the effects of the top-bottom 
weak isodoublet. Indeed, besides the latter QCD corrections and the well-known electromag- 
netic corrections |J, these genuine electroweak effects are by far the most sizeable: top quark 
induced radiative corrections appear in the self-energies of the weak vector bosons || and in 
the coupling of the Z boson to bottom quarks M and lead to contributions which are quadrat- 
ically proportional to the large mass J?]] of this yet unobserved particle. Two types of higher 
order corrections to the previous quantities have been calculated so far. The first one involves 
virtual Higgs boson exchange which leads to potentially large 0{a 2 m^ / M^r) contributions 0. 
Motivated by the large value of the QCD coupling constant, the second type of corrections are 
the two-loop 0(aas) contributions which, for the transverse components of the vector bosons 
self-energies have been calculated in Ref. [9-13] and extensively discussed in Ref. [14-16], and 
for the Z boson coupling to bottom quarks derived recently in Ref. |T7 . 



However, due to the great technical difficulties that one encounters when calculating at this level 
of perturbation theory, all the previously mentioned results involve some simplifying approxi- 
mations. For instance, the Higgs exchange contributions to both the weak bosons self-energies 
and to the Z-bottom quarks vertex as well as the QCD corrections to the latter vertex, were 
calculated in the limit where the top quark is much heavier than the weak vector bosons, an 
approximation which might turn out to be very poor in view of the lower limit on the top quark 
mass [0. A better approximation has been made in the case of the QCD corrections to the 
vector bosons self-energies since results are available in the limit of a vanishing bottom quark 
mass; this is certainly a good enough approximation for practical purposes since the bottom 
quark mass amounts to, at most, a few percent of the top quark mass. 

From a theoretical viewpoint, however, it would be desirable to have a result with no approxi- 
mation involved except, of course, for the one due to the truncation of the perturbative series. 
This would place analytical calculations in the Standard Model at a level comparable to what 
is known in QED where, in a pioneering work, the authors of Ref. [13 have derived exactly the 



vacuum polarization function of the photon at two-loop order. In addition, an exact calculation 
might turn out to be mandatory in the case of a fourth generation, the existence of which is 
still allowed by present experimental data if the associated neutrino is heavy enough ffl. 
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In this article, we present exact and compact analytical expressions for the 0(aas) contribu- 
tions of quark pairs to the vacuum polarization functions of the electroweak gauge bosons in 
the most general case: real and imaginary parts of both the transverse and longitudinal com- 
ponents, for different and non-zero quark masses and for arbitrary momentum transfer. 

The paper is organized as follows. In the next section we summarize all the one-loop results 
which will be relevant to our discussion. In section 3, we discuss the dependence of the result on 
the definition of the quark masses and provide the material necessary to derive the vector bosons 
self-energies in any renormalization scheme. The expressions for the real and imaginary parts 
of the vacuum polarization function in the most general case are given in section 4. In section 
5, we display the expressions of the self-energies in special situations of physical interest: equal 
quark masses, one massless quark and at zero or very high momentum transfer. Finally, section 
6 contains our conclusions. In the Appendix, we will list for completeness the expressions of 
the one and two-loop scalar integrals that we encountered in this calculation. 

2. One— loop results 

To set the notation and for the sake of completeness, we rederive in this section all the one-loop 
results which will be relevant to our discussion. 

The contribution of a fermionic loop to the vacuum polarization tensor of a vector boson i, or 
to the mixing amplitude of two bosons i and j, denoted IP^, is defined as 

0> (2.1) 

where T* is the covariant time ordering product and q the four-momentum transfer; J* , J-j are 
fermionic currents coupled to the vector bosons i,j and constructed with spinor fields whose 
corresponding masses are m a ,m b . The definition of corresponds to +i times the standard 
Feynman amplitude. The vacuum polarization tensor can be decomposed into a transverse and 
a longitudinal part, 

n^ 2 ) = (g»„ - q -f^ iW) + n*V) (2.2) 

and the two components can be directly extracted by contracting IT^, (g 2 ) by the two projectors 
9ixv — CLtiCLv/q 2 and q^qu/q 2 ] in n = 4 — 2e dimensions one has 

i4(g 2 ) = " ^) W) ( 2 - 3 ) 

niV) = ^?n^(g 2 ) (2.4) 



IF 



(« 2 ) 



d 4 xe iq - x 



< 0|T* 



x 
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The one-loop transverse and longitudinal components can be written as [with s = q 2 ] 

OL V .... -i 

U^ L (s) = -N c (v l v 3 + a l a 3 )s L (s) + (uV - a* a 3 ) m a m b lip L (s) (2.5) 

with A^c the number of colors and v l and a 1 the vector and axial-vector couplings of the gauge 
boson i to the fermions expressed in units of the proton charge e = ^/Ana. 

The vector and axial-vector components of the vacuum polarization function [with the coupling 
constants factored out] are then simply given by 

u t£( s ) = sII t,l( s ) ± m a m b U.T jL (s) (2.6) 

which exhibits the fact that U A ' v (s) can be obtained from U v ' A (s) by simply making the 
substitution m a (m b ) — > —m a (—m b ) as expected from 75 reflection symmetry. 

At the one-loop level and for arbitrary fermion masses, m a 7^ m b 7^ 0, the vacuum polarization 
amplitude, Fig. la, reads 



?J ~ C ^ 4tt y J (2tt)" (F-m2)[(A;-g) 2 -m2] 



where /x is the 't Hooft renormalization mass scale introduced to make the coupling constant 
dimensionless in n = 4 — 2e dimensions; we have also introduced an extra term (e 7 /47r) e [7 is 
the Euler constant] to prevent uninteresting combinations of log47r, 7 • • •, in the final result. 
After contracting the amplitude with the tensors {g^ v — q^q u /q 2 ) and q^q^/q 2 , and expressing 
the scalar products of the momenta appearing in the numerator in terms of combinations of the 
polynomials appearing in the denominators, one is led to the calculation of two scalar one-loop 
integrals which, for completeness, are given in Appendix. 

The expressions of 11^ £ (s) in the general case m a 7^ mj 7^ are then 

„. . , 2 + 3a + 3(3 1/2 A , ,5 1, 1 , „ x2 

n?(s) = -- l [^ + a + (3j(p a + p b ) + - + -(a + (3)--(a-(3) 2 

- — (a - (3) 3 log^ + — [3(1 + a + (3) - A] (logx a + loga; 6 ) 
nl(s) = ^-\(a + (3)(p a + p b ) + a + (3 + l(a-(3) 2 + ^(\-l)(a-P)\og^ 
+ ^A5(A - 1 - a - (3)(\ogx a + logx fe ) 

= n £ (s) = ---(p a + p 6 ) + 2 + -(a-/?)log^ + -A^(loga; a + loga; fe ) (2.8) 
where we use the variables [ /3, p& and x b are defined similarly] 



A = 1 + 2a + 2(3 + (a - (3) 2 



(2.10) 



From these general formulae, it is straightforward to derive the expressions of the self-energies 
in two special situations of physical interest: i) the two quarks have equal masses, which cor- 
responds to the case of a neutral gauge boson [the photon or the Z boson] and ii) one of 
the quarks has zero-mass, which corresponds to the case of a charged gauge boson with one 
quark much heavier than the other [as for the top-bottom isodoublet]. For completeness, we 
exhibit in the following the expressions of the vector and axial-vector parts of the transverse 
and longitudinal^ components of the gauge bosons self-energies in these two limiting situations. 

For mj = m a , one has 



n£(*) 



s 
3 
s 
3 


2sa 



4a 



I_, o + 5_ 4£l + (1 + 4a)1 / 2(1 _ 2a)log _ + ^_ )2 



- - Pa ) (1 + 6a) + | + 8a + (1 + 4a) 3/2 log ■ 



--p a + 2 + (l+4a) 1 / 2 log- 



4a 



1 + VT+Iaf 



4a 



1 + ^1 + 40)2, 



(2.11) 



while for m& = 0, one has 

s 



n T («) 



n r . is) 



(2 + 3a) ( - - p a ) + ^ + 2a - a 2 + (2 + 3a - a 3 ) log — — 
\e / 3 1 + a 

-n Pa + 2 + a + (1 + a) 2 log - 

e 1 + a. 



(2.12) 



In many cases, it is of interest to evaluate the vacuum polarization functions at zero momentum 
transfer. This is the case for instance of the radiative correction to the p parameter || which 
measures the relative strength of the charged to neutral current at q 2 = 0. In this limit, the 
contributions of a fermionic loop to the transverse and longitudinal components simplify to 



n$(o) 



ill 

- — (m 2 + m 2 b ) - -{m 2 a + m 2 b ) + -(m 2 + m 2 b )(p a + p b ) 



1 m a + m, 
4 m 2 r , — m? 



mi 



1 

- + 1 

e 



1 1 m 2 + m 2 

2 2 — 



mi 



(2.13) 



Note that in this limit, the longitudinal and transverse components of the vacuum polarization 
function are equal and that for a vector current with equal fermion masses m a = m [as it is 

^^Note that the contribution of the longitudinal part in processes involving fermions in the initial and/or final 
states is always suppressed by the masses of these fermions. 
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the case for the photon] the vacuum polarization function vanishes. 



In the opposite limit m 2 a ,m1 -C s, the vector and axial-vector components of the gauge bosons 
self-energies simply reduce to 



K V L A (s) 



S 

3 




1 , -a 5 
- - log — + - 

e j2 z 6 



(2.14) 



One sees that Hj, L (s) 



Il^ L (s) as expected from chiral symmetry and that, in this limit, there 



is no mass singularity. Note that U?' (s) are quadratically divergent for |g| — > oo. 



In all the previous expressions the momentum transfer has been defined to be in the space-like 
region, s < 0. When continued to the physical region above the threshold for the production 
of two fermions, s > (m a + m&) 2 , the vector boson self-energies acquire imaginary parts, the 
latter being related to the decay widths of the vector bosons into fermions. Adding a small 
imaginary part — ie to the fermion masses squared, the analytical continuation is consistently 
defined. 



From the expressions eqs. (2.8) 
making the substitution 



the imaginary parts can be straightforwardly obtained by 



log x a . 



log \x a h\ + in 



and one has 



Xmlij. (s 



Xmll 



= -A 2 
2 

7T 1 . 
LV*)= o Al t A 



1 + a + P) - -A 



l-a-P] 



ImII T (s) 



XmUr (s) 



7rA2 
71 



(2.15) 



or, equivalently, 



XraU^ A {s) 



XmIlY ,A (s) 



77 ^ 

— SA2 

2 

2 SA2 



;i + a + /3) - -A ± 2 



A-l-ot-^ ± 2 



m a m b 



(2.16) 



The knowledge of the imaginary part of the polarization function, which can be calculated 
directly using Cutkosky rules fL9f , allows an alternative way for obtaining its real part IT(s): 
the latter can be expressed as a dispersive integral of Imll(s). The connection between the 
dispersive approach and the results which are derived using dimensional regularization, has 



been discussed in Ref. |15[ in the special cases m a = rrib and mb = and more recently |16| in 
the general case. 
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3. Renormalization and scheme dependence 

At 0(aas), the two-loop diagrams contributing to the vacuum polarization function EP£, (g 2 ) 
induced by quark loops [up to a factor +i] are shown in Fig. lb. In the 't Hooft-Feynman gauge, 
using the routing of momenta shown in the figure and following the notations introduced in the 
previous section, one can write the bare amplitude as 



bare 



--aa 5 (167r 2 ) 



47T 



2c 



d n k x 



(2vr) n J (2tt) 



d n ko 



(3.1) 



with 



A' 



Tr 



Tr 



(#i + m a )7 M (^ - a t 7 5 )(^ 1 - (/ + m b )j x (^2- 4 + m b )^ u {v 3 - a J 7 5 )(y 2 + ^a)7 ; 

(A;i - /c 2 ) 2 (A; 2 - to 2 ) (A$ - m 2 a ) [{k x - q) 2 - m 2 ] [(k 2 - q) 2 - m|] 
(#i +m a )7„(f i -a i 7 5 )(^i- + m b )'y u (v i -a J 7 5 )(#i + mahxifc + ™a)7 A 



(*1 



k 2 ) 2 {k\ - m 2 a ) 2 ik\ - ml) [(h - q) 2 - m 2 b 



(3.2) 



This bare amplitude has to be supplemented by counterterms. By virtue of the QED-like Ward 
identity, the vertex and fermion wave function counterterms cancel each other and only quark 
mass renormalization has to be included. The latter is obtained by considering the diagram 
shown in Fig. 2a, the amplitude of which reads in dimensional regularization 




d n k 7 A (j/— ^ + m)7A 
(27r) n [(p - k) 2 - m 2 } k 2 



(3.3) 



where p is the four-momentum of the quark and m its bare mass. This expression can be 
decomposed into a piece proportional to (p 1 — m) which will enter the wave function renor- 
malization and another piece proportional to m which will give the mass counterterm. After 
integration over the loop momentum, the latter is given by 



v ( 2\ a s 

^m{P ) = 

7T 



m 2 



m 1 - 2e/3 
7 1 - 2e 



r(i + e) 



1+ 1^-1 )0(e) 
\ m z / 



(3.4) 



The mass counterterm will now depend on the renormalization procedure; i.e. on the definition 
of the quark mass. For instance, in the on-shell scheme which is usually used to calculate 
radiative corrections in the electroweak theory |20f , the fermion masses are defined at p 



m 



and correspond to the position of the pole of the fermion propagators. They are referred to as 
the on-shell or physical masses and the counterterms will be given by 



5m = m(m 2 ) 



m 




1 - 2e/3 
1 -2e 



(3.5) 
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One then inserts this mass counterterm in the one-loop self-energies, as depicted in the dia- 
grams of Fig. 2b, which is equivalent to calculate 



CT 



-8m a 







8m a 



1— loop 



d 



duih 



n^ 2 



1— loop 



(3.6) 



where the one-loop vacuum polarization function is given by eqs. (2.7-2.8). The renormalized 
two-loop self-energies will then read 



re(g 2 ) = re(g 



bare 



+ n^(g 2 ) 



CT 



(3.7) 



However, one can also employ a different definition for the quark mass. For instance one can 
use the MS scheme in which the mass is defined by just picking the divergent term in the 
expression of E m (p 2 ) [the related constants log47r,7, ■ ■ • have already been subtracted] or the 
running mass^. Having at hand the result of the vacuum polarization function in the on-shell 
scheme, one can obtain the polarization function in any other scheme X: one simply has to add 
to the expression of the two-loop self-energy eq. (3.7) the quantity^ 



AIC(g 



x 



Am xmad 



a dm n 



KM 1 



1— loop 



1— loop 



where in terms of the quark mass m^ b defined in the scheme X, one has 



{3.1 



Am y n 



1 - 



rn 



x 

o.6 



(3.9) 



Decomposing the one-loop function into its transverse and longitudinal components as in 
eq. (2.5), one obtains for the derivative of these components 



1— loop 



d 



dm a 



= - - ^(Pa + Pb) + a (P - «) - 2 a ( a ~ lo g ~p 

+ ^ a(l - a + p)(l - 2p\' 1 )(\ogx a + logx fe ) 
a 



1— loop 



+ 



1 la 
- - -a(p a + p b ) + a(2 -3(3 + 3a) + -a[3(a - (3) 2 + 4a] log - 

1,1 
-A 2 a 
2 



1 + 3a - 3(3 + 2(3\~ l (l + (3 - a)] (\ogx a + logx 6 ) 



2 A discussion of the phenomenological consequences of using a different definition for the quark masses will 
appear elsewhere pl[ . 

3 Note that since at this stage we are only discussing the difference between two renormalization schemes, we 
do not need the 0(e) terms in the one-loop result of the vector bosons self-energies. These terms are of course 
needed to evaluate the renormalized two-loop self-energies in a given scheme. 
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m«7 — m a m b ]l T L (s, 



1— loop 



= m Q m 6 
+ -A. 



11 1 (y. 

- - ^(Pa + p b ) + 2 + -{3a- (3) log ^ 

1 + 2aA~ 1 (l + a - 0) (log x a + log x b 



(3.10) 



and similarly for the piece involving m& which can be obtained by making the substitution 
m a <-> mj,. Hence, once the complete result of the vacuum polarization function is known in 
the on-shell scheme, it is straightforward using eqs. (3.8)-(3.10), to obtain the corresponding 
results in any renormalization scheme. At this stage, a few remarks are mandatory. 

First, one notices in the previous expressions the occurence of terms that are inversely propor- 
tional to the velocity factor A 1//2 . In principle, after an analytical continuation to the physical 
region beyond the threshold for quark pair production, s > (m a + m&) 2 , these terms would 
diverge for energy values near the production threshold, A ~ 0. However, as we will see later, 
when evaluated in the on-shell scheme, which is the only scheme where the physical threshold is 
well defined, the renormalized two-loop polarization function is free of these A~ 1//2 factors. Near 
threshold, the dominant terms will be constants and would correspond, once the vector boson 
self-energy is normalized to its one-loop value, to the well known A -1 / 2 Coulomb singularities 
which require a non-pertubative treatment; see for example Ref. ||15|| . Of course, the A -1 / 2 
terms can be present in the vacuum polarization function when it is evaluated in a different 
scheme but in this case, the threshold is not well defined since the masses are not "physical" 
masses. 

In principle one can define the masses m a and m b in two completely different schemes; this 
can be useful if, for instance, different scales are involved in the evaluation of the vacuum 
polarization functions. Taking the example of the top-bottom isodoublet, one can employ the 
on-shell mass for the heavy top quark, which is suitable if one wants to discuss threshold and 



possibly non-pertubative effects [ 15| ; and use a running mass evaluated at the scale q for the 



relatively light b quark, which in general, avoids the appearence of large logarithms for q 3> ml 
||22|| . In the rest of our discussion, however, we will stick to the on-shell mass scheme. 

Finally, the scale at which the strong coupling constant as is evaluated can, of course, be 
completely different from the one chosen for the quark masses. In general, as is evaluated at 
the scale of the problem at hand, i.e. as = as{s). However, for heavy virtual quarks, the 



effective coupling is usually taken at the mass of the quark [23] and this choice can be justified 
by arguments based on effective field theory p4j [for a discussion on this point see also the 

0] 



second paper of Ref. 



For instance, in the case of the top-bottom isodoublet, the two-loop 



QCD corrections to electroweak observables are calculated with as(m 2 ) while for light quarks, 
one usesQ as(q 2 ). 



An exception is the contribution of the light quarks to the photon self-energy which often has to be evaluated 
at zero momentum transfer. Since both q 2 and m 2 are small compared to the QCD scale, one can no more use 
perturbation theory. However, the problem can be circumvented by using a dispersive approach and relating 
the photon self-energy to e + e~ — ► hadrons low energy data plfl. 
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4. Exact two— loop results 



In this section we give the expression of the vacuum polarization function at order 0(aas) in 
the general case m a 7^ m;, 7^ and for arbitrary momentum transfer. The result will be given in 
the on-shell mass scheme. We will follow very closely the notations and definitions introduced 
in section 2; since confusion should be rare, we will use the same notation for the one-loop and 
two-loop self-energies and m a , b will stand for the on-shell masses. 

At 0(aas), the transverse and longitudinal components of the vacuum polarization function 
Yl^ u (q 2 ) are defined by [the color factor Nc is now included] 



m L {s) 



a as 

TT TT 



(uV + aV> n+ L ( s ) + 



l - J -^m a m b 11? As) 



[V V 



a a 



(4.1) 



where IT^ L are the sum of the corresponding components in the bare two-loop amplitude 
eqs. (3.1-3.2) and the mass counterterm which can be obtained from eq. (3.11). 

Similarly to the one-loop case, after contracting EP£, in eqs. (3.1) by the tensors {g^ v — q^q 11 / q 2 ) 
and q^q v jq 2 and expressing the scalar products of momenta appearing in the numerators in 
terms of combinations of the polynomials in the denominators, one is led to the calculation of a 
set of scalar two-loop integrals. Most of these integrals have been first calculated by Broadhurst 
in Ref. f26j; the remaining integrals reduce after straightforward computations to the previous 



ones |27], [TO] . For completeness, these scalar integrals will be given in the Appendix. 



After a cumbersome computation, and taking advantage of the symmetry in the change a. *->■ 3, 
we can write U? L (s) in a relatively simple and compact form 



Ih1 



3 " 1 f l o 11 \ !/ W11 1 n ON 55 71 

~2^ + 7 V4 + ~ ~4 a ) + 4 + ~ ~ a + + 24 ~ 24 a 

112 1 

-a 2 H a3 + -a[G(x b ) - G(x a )} + - log x a log x b [3(a + 3) + 2a3(A + a + 3)] 



6 



6 





1 


+ 




12 




1 


+ 




12 






+ 


V 





- \ogx a [(a -13+ A2) (11 + 19a + 19/3 + 12a/3 - 5A) + (a - 3)(42 - 5a - 53) 
log 2 x a [(1 - 3a - 38) (A - 1 - a - 3 + (a - /3)A^) - 9(a + 3) + 8a/3 



7T 



[{a + 3 — 2)A - 12a/3] I -- [3{1 + a + 3) - A] T 

3 



+ {a< >3} 



n 



T 



2e 2 



11 



4 



11 



11 



7T 



77^ + - 3p a - — + —p a - -(p a + p b ) + — + a - — 



4 



+ ilogx a [a-3+(a + 3 + 9)(a-3 + X^ 



log 2 x a 1 + a + 3 - A + (3 - a) A 



+ - logx a \ogx b [3(1 + a + 3) + 2a3\ -{\ + a + 3)Z- 2J'| 



+ {a< — >3} 
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T-T+ f 3a 1 / 11 \ a . w , 3 7 , 13 „ 

HI = + - (^3ap a - — aj + -(p a + p fe ) (11 - 9p a + 3p 6 ) + -a + -a - ya/? 

(a - /3 + A*) (x-l-a - P- ^ a P^j + {a - (3) (a + (3 + 



+ - log x a 



20 
~9~ 



7T 



+ - log 2 x a [A(2A - 2 - a - /3) - 3(a + p) - 16af3 + (2A + 1 + a + (3){a - f3)\^ 
+ - log x a log x b {a + /3)(1- 2a/3) + ^ [4a/3 - (a + P)X] l+{l + a + (3- A)J'| 



+ {a 







n; 



[3 1 / 11\ 11 3 , , 2 3 7T 2 

3 i 1 i 

- -logx a a-/3 + (a + /3-3)(a-/3 + A2) + -log 2 x a l + a + /3-A + (/3- a)A^ 



3 1 
+ - log x a \og x b [l + a + (3 -2af3]-{l + a + (3)1 -2T\ + {a 



0} 



(4.2) 



with 1 and X' given by 

1 = F(l) + F{x a x b ) - F{x a ) - F{x b ) 

T = \^G{x a x b )-^{(3-a + \i)G{x a )-^{a-f3 + \^)G{x b ) 



(4.3) 



In terms of the polylogarithmic functions ^(x) = — J y 1 log(l — xy)dy and Lia(x) 
— Jq 1 y^ 1 logy log (1 — xy)dy, the functions F and G are given by 

- 6Li 3 (x) — 4Li 2 (a;) logx — log 2 xlog(l — x) 



F(x) 
G{x) 



2Li 2 (x) + 21ogxlog(l — x) + 



x 



1 — X 



log X 



These two functions also admit a simple and useful integral representation 



F(x) 



\ogy\ x 
dy log - 

" \ l -yJ y 



G(x) = x F\x) = [ X dy( 

Jo \i-y 



(4.4) 



(4.5) 



The imaginary parts of the vacuum polarization function are derived along the same lines as 
discussed previously in the one-loop case. Using the fact that 



IvTL log X a h = 7T 



1ml' 



x 



1m log 2 x atb = 2n log \x a , b \ 

-(/3-a + A5) [log(l-x a ) 



log(l - x a x b ) + XaXb log \x a x b \ 



1 — X, 



logical 



1 X a %\) 

log(l - x b ) + 



x b 
1 - x b 



\og\x b \ 
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XmX = ix J = -2n [ALi 2 (x a x b ) - 2Li 2 (a; a ) - 2Li 2 (x fe ) + 21og|x a Xfe| 
x log(l - x a x b ) - log \x a \ log(l - x a ) - log \x b \ log(l - x h )\ 



(4.6) 



one obtains for XrnlY^ L 



—XmWp{s) 
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XmH T (s) 



+ 



~(11 + 19a + 19/3 + 12a/3 - 5A)A^ + -(a - 0) [log(l - x b ) - log(l - x a )} 
6 3 

3(1 + a + /3)((/3 - a)A5 + 1 + a + /3 + 2a/3 - A) + 8a + 26a/?l log |ar | 



2 
3 



(9 + a + /3)A* + [4(1 + a + p) + 2a/3 — X + (J3 — a)\*\ log \x a \ 
4(1 + « + /?) + 2a/3 - A + (a - /3)A*1 log |x 6 | - 2(1 + a + /3) J" - 4J"' 







+ 




6 t 




1 r 


+ 




6 V 






+ 





7T 



Jmn+(s) 



+ 
+ 



(9A - 9 - 9a - 9/3 - 12a/3)A2 - [(a + /3)A - 4a/3]^ + 2[1 + a + (3 - \\J' 



2 

1 r 



(2A+ l + a + /3)A5(a-/3 + A2) - (A + 3a/?)(2a + 2/3 + 3) - 7a/3 log|x a 
(2A + l + a + /3)A5(/3-a + A^) - (A + 3a/3)(2a + 2/3 + 3) - 7a/?] log|x fe 



7T 



-XmIl L (s) 



+ 



3(3 - a - /3)A5 + [4(1 + a + /3) - 6a/3 - A + (/? - a)A*] log |x a | 
4(1 + a + /?) - 6a/3 - A + (a - /3)A*1 log |x 6 | - 2(1 + a + /3) J - 4,7' (4.7) 



These expressions for the real and imaginary parts are the main result of this paper. 



Close forms for XmU^' L (s) in the general case m a ^ m b ^ have been also derived in the past 



by a number of authors [IT], |29|, [3(| [in the first reference only the transverse part is given] by 
directly calculating the QCD corrections to the flavor changing decay of a vector boson. The 
results that we obtain here using a completely different method, agree with those of Ref. JTlJ 
and also with Ref. [B0J once some obvious mistakes in the integrals of their Appendix [Ji and 



J 2 ] are corrected; see also Ref. |31] 
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5. Special Cases 

From the general formulae given in the previous section one can derive the expressions of the 
vector and axial-vector components of the self-energies in various special situations of physical 
relevance. Here we will exhibit these expressions in four different cases: i) the two quarks have 
equal masses, %%) one of the quarks is masslessQ and Hi) the quark masses are much larger or 
much smaller than the momentum transfer. The results that we obtain in these special cases 
provide several checks of the general expressions of the previous section and allow for a com- 
parison with various results available in the literature. 



i) Case nib = m a 

The real parts are given by 



l 55 26 2 

M p H a + \/T+ 4a (1 — 6a) log a; a(4 + a) log 2 x 



2e 



12 



+ -(4a 2 -1) \f{1)+F{x 2 )-2F{x) 



[1 - 2a) VI + 4a \G(x 2 ) - G(x 



{6a 1 55 19 
+ (1 + 24ap - 22a) (1 + 12ap - 22a)p + a + 4a 2 - avr 2 
e 2 2e 12 6 

2 2 

+ (1 + 12a + 4a 2 ) Vl + 4a log x + -a(5 + 11a + 6a 2 ) log 2 x (1 + 2a) 

3 3 

x (1 + 4a) \F(1) + F(x 2 ) - 2F(x)} - ^(1 + 4a) 3/2 \G(x 2 ) - G(x) 



f 3 / 11\ 1 3 1 9 / x / , 

2sa|--+ (6p-yJ - + - - 6a - -7T 2 + 3(3 - 2a) Vl + 4a log x 

+ (11 - 6p)p + (3 + 4a - 6a 2 ) log 2 x - 2(1 + 2a) [F(l) + F(x 2 ) - 2F(a 
- 4V1 + 4a [G(x 2 ) - G(x)] } (5.1) 
with x = 4a/(l + Vl + 4a) 2 , and the imaginary parts are 

ys j~A^(l - 6a) - 2a(4 + a) log |x| - (1 - 4a 2 ) J - 2(1 - 2a)J'\ 
y s [^H 1 + 12 « + 4a 2 ) + 2a(5 + 11a + 6a 2 ) log |x| - (1 + 2a)AJ - 2\J'\ 
4vrsa|^A5(3-2a) + (3 + 4a - 6a 2 ) log \x\ - (l + 2a)J-2J'| (5.2) 



Xmlly (s 



where 



J 
J' 



Li 2 (:c 2 ) — 2Li 2 (x) + 2 log |x| log(l — x 2 ) — log \x\ log(l — x) 



2V1 +4a 



2 logf 1 - x 2 



, /■, \ x (3x — 1) , , , 

log(l — x) H 5— log x 

1 — x z 



(5.3) 



5 In fact, we have expanded the expressions of the previous sections around mi, = m a and nib = retaining 
terms up to order (mb — m a ) 2 /s and rn^/s, respectively. The rather lengthy expressions can be found in Ref. 32 
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The vector part of the longitudinal component vanishes in this case: this is expected to occur 
as a consequence of the QED-like Ward identity. 

In the equal mass case, Ily(s) coincides, up to a color factor, with the irreducible part of the 
photon self-energy in QED which has been calculated in the fifties by the pioneers of Ref. ]|il| . 
Ut' A (s) has also been derived using a dispersive approach in [II], and more recently in [13|| . 

ii) Case = 

In this case, the coefficients of II^ T (s) vanish so that II^ r (s) = nf T (s) = sll^ r (s), and one 
obtains for the real parts 



EV (s) 



11 



1 



2e 



11 



55 71 



12 24 



+ 



-of + (1 



1 3 
l + -a 



\ogx + -(1 




3a) (1 + a) log 2 a; 
^(2 +a)(a- l)G{x) + i(a - 2)(1 + a) 2 [F(l) - F(x)]} 



sa 



1 



2e 2 



11 



4 



7 



11 



7T 



9 



— ^ + — - + 3p - + - + -a + — p - 3p 2 - — + -(1 + a) 2 logx 



+ ^(3 + 2a)(l + a) 2 log 2 x + (l + a)G(x 
with x = a/(l + a), and for the imaginary parts 



)-(l + a) 2 [F(l)-F(x)]} (5.4) 



XmII^' A (s) 



TmUY' A (s) 



( 5 2 5 1 a 

7ts U + -a + -a 2 a 3 (1 + a) 2 (4 - 5a) log(l + a) (5a 2 + 4a - 5) 

1 2 3 6 3 3 



x log(-a) 



;(2-a)(l + a) 



2Lio 



a 



a + 1 



- log(l + a) log 



-a 



-7rsa(l + a) 



+ 2(1 + a) 



2Li, 



1 + a) 

a 
a + 1 



9 

-- + (5 + 2a) log(l + a) 



1 + aj 
(3 + 7a + 2a 2 ) log(-a) 



log(l + a) log 



-a 



1 + a 



(5.5) 



We see that in this limit, the self-energies are free of mass singularities as required by the 



Kinoshita-Lee-Nauenberg theorem [03]| . The expression of Ely (s) in this special case has 
been obtained in Ref. [13] by means of a dispersion integration. 

Note that the results for the real parts of the vector bosons vacuum polarization function in 
these two special cases slightly differ from those obtained some time ago by one of the present 
authors in |1(J: the it 2 terms are absent in the expressions given in the latter reference [this is 
due to the fact that the (l+e 2 7r 2 /12) factor in the mass counterterm eq. (3.5) has been omitted^]. 
However, this discrepancy does not affect the radiative corrections to physical quantities when 



°We thank B. Kniehl for bringing this point to our attention. 



14 



evaluated in the on-shell scheme, which is commonly used to calculate radiative corrections 
to electroweak observables. Indeed, these ir 2 terms are related to the 1/e 2 poles which must 
cancel in physical quantities and therefore they do not appear in the radiative corrections to 
observables like the p parameter, the W boson mass or the effective weak mixing angle sin 2 9w 



Hi) Cases s 3> m\ b and s <C m 2 



b a,b 



In the limit of zero momentum transfer, the expressions of n^'y(s) simplify considerably, and 
one obtains [we use the fact that G(x) + G(l/x) = 27r 2 /3] 



n?2(o) 



3 . o 2\ 

— (m a + m b ) + - 



11 



m 2 + m 2 b ) - 3m 2 a p a - 3m 2 p fe 



11 



rnlpa + m 2 b Pb) 



+ 3m 2 p 2 + 3m^ 2 + — (m 2 + m 2 ) + -(m 2 



m. 



G\ — 



\ m 2 , 



+ (m 2 +m 2 )^+ , 
4 mi 



mlml , m 2 

' log —I + m 2 m 
m£ m£ ( /// 



2 2 m a + m b i 2 
2 ,,,2 / ...2 _ ,,.2\2 l0 § „,,2 



m b) ,, % 

- \3p a + 3p b - —j - -{p a + p b ) + 4(p a + Pb ) - — 



+ 3 ™fe ~ m \ Pa + 3- 



m 2 m 2 



m 



"6 \ m a ~ "%) "°b * 

As expected, the self-energies are regular in this limit and one has again IT^' A (0) = 11^(0). 
In the two special cases m a = m b and m& = this expression simplifies to [HV(0) vanishes in 
the equal mass case] 



m 



2\2 



log 2 



mi 



1 



7T 



(5.6) 



n^ L (0,m a = m b ) 



m„ 



n^(o,m fe = o) = -f 



6 



6 



1 31 
-(11 - 12p a ) + 12p 2 - 22p a + — + 7T 5 

C — 



1 



35 



- + -(11 - 12p a ) + 12p 2 - 22p a + — + -7T 



(5.7) 



Note also the presence of the extra 7r 2 terms compared to Ref. |TtJ as discussed previously. 
In the opposite limit \s\ ^> m 2 fe , the vector and axial-vector components are simply 



~ s 
~ 



1 

27 



log — + 

p 2 



55 
12 



4C(3) 



(5.8) 

with £(3) = F(l)/6 = 1.202, in agreement with Ref. J|34|] . Here again, there is no mass singu- 
larity and the quadratically divergent components I1 T ' (s) obey the chiral symmetry relation. 
The next term of this expansion, i.e. the terms of 0(m 2 b /s), can be found in Ref. [|16 |. 



In view of the several checks that we have performed and of the various results available in the 
literature that we recover in the special situations discussed above, we have good confidence 
that the expressions that we give here for the general case m a ^ m b ^ are correct. 
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6. Summary 



In this paper, the contribution of heavy quarks to the vacuum polarization function of the 
electroweak gauge bosons, and to the mixing amplitude of two gauge bosons, were calculated 
at first order in the strong interaction. Both the transverse and longitudinal components of the 
vacuum polarizations functions were given in the most general case: arbitrary masses for the 
quarks and finite momentum transfer. Full analytical formulae for the real as well as for the 
imaginary parts at O(aas) were presented in the on-shell quark mass scheme. They are given 
in a compact form in terms of polylogarithmic functions which can be easily evaluated using 
standard computer routines. 

The dependence of the result on the definition of the quark masses has been discussed in detail. 
We provided the necessary expressions which allow, having at hand the result in the on-shell 
mass scheme, to obtain the vector bosons self-energies for quark masses defined in any other 
renormalization scheme. 

The expressions of the self-energies were then derived in two situations of great phenomeno- 
logical interest: first the case where the two quarks have equal masses which corresponds to 
the photon and the Z boson self-energies and second, the case where one of the quarks has a 
negligible mass with respect to the other which corresponds to the approximate contribution 
of the top-bottom isodoublet to the charged W boson self-energy. The expressions of the self- 
energies for asymptotic values of the momentum transfer: s — > oo or s — > were also given. 

Some of the results in the special cases that we have considered here for the real parts of the 
vacuum polarization functions, as well as the imaginary parts in the general case, are available 
in the literature and the comparison with them provided very powerful checks of our calculation 
in the general case. 
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APPENDIX: Scalar Integrals 



In this appendix we will give for completeness the expressions of the one and two-loop scalar 
integrals which are needed in the evaluation of the two-loop vacuum polarization functions 
26], [27|, 1Q| . We will closely follow the notation of Ref. |Hj], but simplify further some of the 



expressions and correct a few misprints which occured in the Appendix. 
One— loop integrals. 

These are needed in the calculation of the amplitude both at the one-loop and the two-loop 
levels. Indeed, when expressing the scalar products of momenta appearing in the numerator of 
eq. (3.1) in terms of combinations of the polynomials in the denominators, the gluon propagator 
[i.e. the {k\ — k 2 ) 2 factor] cancel in some cases and the two-loop amplitude reduce to the product 
of two one-loop integrals. These one-loop intregrals are 



d n k 



71 



i/2 Ul 



mi 



-a 



1 1 1 A, 2 ! 1 ^ 

- + 1 — log a + e I - log a - log a + 1 + — 



(a.i; 



[and similarly for Db where one has to replace a = —m 2 a /q 2 by (3 = —m 2 /q 2 ] and using the 
variables x a ^, A defined in eqs. (2.9) and (2.10) 



K 



-q 2 e^Y 



d n k 



Tj-n/2 ^2 _ 7772][(& — q) 2 — 77t|] 



- + 2 + (a - (3 + Aa ) log x a - log a 



- [a <- 0\ + 0(e) 



(A.2) 



Note that in the expression of K, the 0(e) term is not needed: its total contribution, via terms 
likes K 2 , KDa-i KdK/da, ■ • • , vanishes in the final result. 

The derivatives of K and Da with respect to a will enter the expressions of the two-loop 
integrals and are given by 



ODa 
da 
OK 



da 



2\2 



1 — e log a + e 2 f - log 2 a H 

& 1 2 & 12 



[1 + a- (3 + Aa) logx a + (1 + a - f3 - As) logx b 



(A.3) 



Two-loop integrals. 

To simplify the notation, we shall use for the definition of the two-loop integrals the following 
abreviations: 



V ' J 7T™/ 2 J W 2 



(A.4) 
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and 



K h2 = k\ 2 -ml , K = (h-k 2 ) 2 , Q 1>2 = (fc li2 - g) 2 - m 2 (A.5) 
With the help of the variables defined in eqs. (2.9)-(2.10), the two-loop integrals are 

q 2 P = < 

W A = i < ^T- >= 



Ma 
Va 
-Ra 
iV A 

Qa 
Sa 



Ki K 2 ° Qi q 2 >=(!"« + + (! + «- - ^(£>a - ^) 2 - ^ 2 



(g 2 ) K 2 ^ 2 * 9a 

1 0i o , (9-D 4 

1 n <9fT 



1 

r < 



l 



g 2 K^Ko 
1 



>= 



= < 



K 2 K 2 K 



> 



da 
1 1 -e 
~a 1 - 2e 
1 1 - e 



2a 2 1 - 2e 



^4 -D, 



2a- 



dD f 



1 < (fc 1 -2fc a) - g 



g 2 K 1 K 2 K Q 1 
1 



log^q 
1 3T n 



da 



x b \ogx h 



= < 



= < 



wf„o7 >= I (1 + £)A ' 2 + ^ 



3- 



1+ , 

1 — X b 

Xq lOg 2 Xq _ X b lOg 2 X b 

{i-xqf {\-x b y 



- G{x a ) + G{x b ) 



K!Q 2 K 



>- 



I + -(2 + e)K 2 + M A + N A -(l-a + (3)J A 
o 4 



1 



> ^wkir >= Iki-* + P)Ja-l + m a ] + n a 

q z KiK 2 Q 2 K 2 



< 



< K 2 K 2 Q X K« > \i 



lOg Xq 
1 - Xq 



l-X a . 



x b \ogx b ( ^ x b \ogx b 



1 - x b 



1 - x b 



1 

X-K^r- H 

aa 2a 

.2 



A^G(avr b ) - i(l - a + (3 + X^)G(x a ) + i(l - a + /3 - A5)G(a; fe ) 



< 



5 



>= F(l) + F(:r a :r b ) - F{x a ) - F(x b ) 



K x K 2 Q x Q 2 Kq 

The expressions of the functions F and G are given in eqs. (4.4) and (4.5). 



(A.6) 



All the other two-loop integrals can be straightforwardly obtained by taking advantage of the 
symmetry in the interchange of the integration variables k\ and k 2 as well as in the interchange 
of the masses m a and m b . 
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Figure Captions 

Fig. 1 Feynman diagrams for the contribution of quark pairs to the vacuum polarization func- 
tion of a gauge boson at the one-loop (a) and two-loop (b) levels. 

Fig. 2 Feynman diagrams for the one-loop quark self-energy diagram (a) and for the contri- 
bution of the mass counterterm to the polarization function at the two-loop level (b). 
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